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Abstract. In this paper, we consider the problem of the estimation of a 
Weibull tail-coefficient 0. In particular, we propose a regression model, 
from which we derive a bias-reduced estimator of 0. This estimator is based 
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also established. A small simulation study is provided in order to prove its 
efficiency. 
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1 Introduction 



Let Xi, ...,X„hea sequence of independent and identically distributed 
random variables with distribution function F, and let Xi „ < ... < X„ „ 
denote the order statistics associated to this sample. 

In the present paper, we address the problem of estimating the 
WeibuU tail-coefficient > defined as 



1 - Fix) = exp(-H(x)) with H-\x) := mf{t : H(0 > x} = x^£{x\ (1) 



13 investigated this estimation problem and proposed the following 
estimator of 0: 



where A:„ is an intermediate sequence, i.e. a sequence such that 
kn ^ oo and fc„/n — > as n — > cx). 

We refer to lfT5l[T3i l4l and [7] for other propositions and to |^ for 
Local Asymptotic Normality (LAN) results. Estimator ^ is closed 
in spirit to the Hill estimator |[T7| in the case of Pareto-type distribu- 
tions. In [15J , the asymptotic normality of 0„(fc„) is established under 
suitable assumptions. To prove such a result, a second-order condi- 
tion is required in order to specify the bias-term. This assumption 
can be expressed in terms of the slowly varying function £ as follows: 

Ass\xmpiion.{Rc{h, p)) There exists a constant p < Q and a rate function 
b satisfying b{x) — > as x — > oo, such that for all e > and 1 < A < oo, 
we have 



where ^ is a slowly varying function at infinity satisfying 
^^l,asx^oo,forallA>0. 




(2) 



n-i+l,n 

- logX„_|c„+i,„) 



(3) 




sup 

Ae[l,A] 



log{i{Ax)/£{x)) 
b{x)Kp{A) 



1 < £, for X sufficiently large, 
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It can be shown that necessarily \b\ is regularly varying with index 
p [[T4|. Moreover, we focus on the case where the convergence ^ is 

slow, and thus when the bias term in d„{k„) is large. This situation is 
described by the following assumption: 

oo as X oo. (4) 

Let us note that this condition implies p > -1. Gamma and Gaussian 
distributions fulfill dU), whereas WeibuU distributions do not (see 
Tabled]) since, in this case, the bias term vanishes. 

Using this framework, we will establish rigorously in Section|2]the 
following approximation for the log-spacings of upper order statis- 
tics: 

Zj := i log(n/;)( log X„_y+i,„ - log X„_y,„) 



. (e + KH(«ft.))(g^) )/,. (5) 

for 1 < y < kn, where (/i, —,fk„) is a vector of independent and stan- 
dard exponentially distributed random variables. This exponential 
regression model is similar to the ones proposed by ||5l|6| and [[lT| in 

the case of Pareto-type distributions. Ignoring b{\og{n/k„))(^ iog{n/k) ) 
in (O leads to the maximum likelihood estimator 



1 ^" 



" 7=1 



which turns out to be an alternative estimator of The full 

model © allows us to generate bias-corrected estimates d„{k,j) for 6 
through maximum likelihood estimation of 6, fc(logn/fc„) and p for 
each 1 < kn < n - 1. An alternative to this approach consists in using 
a canonical choice for p and to estimate the two other parameters by 
a least-squares method (LS). For the canonical choice of p, we can use 
for instance the value -1, which is the same as the one proposed by 
flTI for the regression model in the case of Pareto-type distributions. 
The asymptotic normality of the resulting LS-estimator is established 
in Section m An adaptive selection method for kn in dn{k„) is also 
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derived. In order to illustrate the usefulness of these results, we 
provide a simulation study in Section H] as well as an application to a 
real data set in Section |5l The proofs of our results are postponed to 
Section [71 



2 Exponential regression model 

In this section, we formalize (O. First, remark that 

F-Hx) = [- log(l - x)ri{- log(l - X)). 

Since X„_y+i^„ = f"^(Lr„_y+i^„), 1 < j < k„, where Uj^„ denotes the /-th 
order statistic of a uniform sample of size n, we have 

Xn-j+l,n = [-log(l - U„-j+Xn)] ^ (- log(l - U„-j+i,n)) 

which implies that 

logX„_y+i,„ = 01og[-log(l - Un-j+l,„)\ + log[^(-log(l - Un-j+l,„))\ . 

Moreover, considering the order statistics from an independent stan- 
dard exponential sample, En-j+i,n = - log(l - Un-j+i,n)- Therefore 

log Xn-j+\,n = 01og(E„_,+i,„) + log[AE„_,+i,„)] 
=: A„{j) + Bn{i). 

Recall that Zj = ; log(n/;) (logX„_y+i,„ - logX„_y,„), 1 < j < K. Then, 
our basic result now reads as follows. 

Theorem 1 Suppose © holds together with {Re{b, p)) and (S]). Then, if 
k„ ^ oo and log k„/ log n — > 0, we have 



sup 

i<i<k„ 



Zj-\d + b{log{n/k„))\ 



login/ j) ^ 
login /kn)f 



Op(&(log(nA„))), (6) 



where ifi,-;fk„) is a vector of independent and standard exponentially 
distributed random variables. 



The proof of this theorem is based on the following two lemmas: 
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Lemma 1 Suppose ([I]) holds together with {Re{b, p)) and (IH). Then, if 
k„ — > oo and k„/n — > 0, we have 

sup |;log(n/;)[A„(;)-A„(; + l)]-0/^| =Op(&(log(n/?c„))), 

i</<fc„ 

and 

Lemma 2 Suppose © /zoZrfs together with {Re{b, p)). Then, ifk„ — > oo and 
log kr,/ log n — > 0, we /zflue 

;log(n/;) [B„(;) - B„(; + 1)] - fc(log(n/fc„)) f 



sup 

'^<i<k„ 



Jog(n/fc„)^ 



= Op (&(log(n/fc„))) . 
The proof of these lemmas is postponed to Section [71 
Corollary 1 Under the assumptions of Theorem^ we also have 



sup 

i<j<k„ 



Zj-\e + h{\og{nlkn))\ 



iog(^/;) 

login /kn)} 



fi 



Op (fc(log(n/;c„))) , 



where {fi,—,fk„) is a vector of independent and standard exponentially 
distributed random variables. 

This implies that one can plug the canonical choice p = -1 in the 
regression model © without perturbing the approximation. From 
model © we can easily deduce the asymptotic normality of the 
estimator 0„(fc„), given in the next theorem: 

Theorem 2 Suppose © holds together with {Re{b, p)) and ©. Then, if 
k„ oo, Vfc^&(log(n/fc„)) A e Kand, if A = 0, logfc„/logn 0, we 
have 



The Asymptotic Mean Squared Error (AMSE) associated to d„{k„) is 
thus given by: 



AMSE{e„{k„)) = ^ + 



fc(log(nA„))^£ 



k„ I 1 / ^p^ 



log(n/;) 



k„-L-J\log{n/kn)l 

1-^ J 



(7) 



This model (O now plays the central role in the remainder of this 
paper. First, it allows us to generate bias-corrected estimates of 6. 
Second, it leads to the number of upper order statistics k„ to be used 
in d„{k„) by minimizing the AMSE given by © after replacing 6, 
b and p by estimators. These two points are described in the next 
section. 



3 Bias-reduced estimates of 6 and adaptive selection 

of kn 

In order to reduce the bias of the estimator OniK), we can either 
estimate simultaneously 0, fc(logn/fc„) and p by a maximum likeli- 
hood method or estimate 6 and by a least-squares approach after 
substituting a canonical choice for p. In fact, this second-order pa- 
rameter is difficult to estimate in practice and we can easily check by 
simulations that fixing its value does not much influence the result. 
This problem has already been discussed in ISHU and |[TT| where 
similar observations have been made in the case of Pareto-type dis- 
tributions. The canonical choice p = -1 is often used although other 
choices could be motivated performing a model selection. 

In all the sequel, we will estimate 6 and fc(log(n/fc„)) by a LS- 
method after substituting p with the value -1. In that case, we find 
the following LS-estimators: 

'dniK) = Zk„ - b{log{n/k„))xk„ 
Z^'UiXj-XkJZj 



&(log(n/fc„)) 



where =(5^)" , x,„ = ^ T!;L, and Z,, = E^l, Z, Our next 
goal is to establish, under suitable assumptions, the asymptotic nor- 
mality of d„{kn). This is done in the following theorem. 

Theorem 3 Suppose ([I]) holds together with {Re{b,p)) and (SI). Then, if 
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kfi ^ oo such that 



login /kn) 



&(log(n/fc„)) ^ A G R and, if A = 0, (8) 



' ana - — n— - oo, (9) 



we have 



login /kn) login /k„) 



login I kn) 



[enikn)-e)^NiQ,e\ 



Remark that the rate of convergence of d„ik„) is the same as the one 

of d„ik„) in the cases where both A and A are not equal to 0. The proof 
of this theorem is postponed to Section [71 

We can also take benefit of the estimation of &(logn/fc„) by esti- 
mating the AMSE given in (O by: 



AMSEiOnikn)) = + 



Klog(n/U)^i:(^^^^ 



Then, the intermediate sequence kn can be selected by minimizing 
the previous quantity: 

kn = arg min AMSE(0„(fc„)). 

This adaptive procedure for selecting the number of upper order 
statistics is in the same spirit as the one used by l|T9l in the context of 
the extreme value index estimation. 

In order to illustrate the usefulness of the bias reduction and of 
the selection procedure, we provide a simulation study in the next 
section. 



4 A simulation study 

First, the finite sample performances of the estimators dnikn), dnikn) 
and dnikn) are investigated on 6 different distributions: r(0.25, 1), 
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r(4,l), |W|(0,1), ^(0.25,0.25), ^(4,4) and D(l,0.5), see the ap- 
pendix for the definition of the latter distribution. We limit ourselves 



to these three estimators, since it is shown in | |15| that dnikn) gives 
better results than the other approaches ['$,^. In each case, N = 100 
samples iP(n,i)i=i,...,N size n = 500 were simulated. On each sample 
{?(n,i), the estimates d„^i{k„), d„j(k„) and d„j{k„) were computed for 
kn = 2,. .., 360. Finally, the Hill-type plots were built by drawing the 
points 



]V ^n,i{kn) I kn, — dn,i{kn) ^nd k„, — dn,i{kn) 



! = 1 



!=1 



( = 1 



We also present the associated MSE plots obtained by plotting the 
points 



!=1 

1 ^ 2 



/ V 



and 



! = 1 



The results are presented on figures HHH In all the plots, the graphs 
associated to d„{k„) and d„{k„) are similar, with a slightly better be- 
havior of d„{k„). The bias corrected estimator 9n{k„) always yields a 
smaller bias than the two previous ones leading to better results for 
Gamma, Gaussian and !D distributions (figures [U-Sl), even though a 
wrong value of p is used (figure SJ. On Weibull distributions, where 
the bias function is zero, (figures |5]-[6l), it presents a larger variance. 
Second, we investigate the behavior of the adaptive procedure for 
selecting the number of upper order statistics in dn{k„). For i = 
1, . . . , N, we denote by 

fc„; = arg min AMSE(dni(kn)) 

k„e[l,350] 

the value selected on the sample {?(n,i)- Note that, as in fT9], in 
our simulations, we limited the range from which kn is selected to 
{!,..., 350}. The mean and the standard deviation of this estimation 
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on the N samples are given by 



N ^ N 



/'=1 1 (=1 

As a comparison, we introduce the value that would be obtained by 
minimizing the true AMSE: 

kT = arg min AMSEidniK)). 

fc„6[l,350] 

On each sample {P(n,i), the estimation of 6 obtained with the selected 

parameter k„^i is given by On,iik„,i). The associated empirical mean 
and standard deviation are: 



N , N 



fi(0,O = QnAKi) and a{dn) = A ]^ (^"'^(^"''^ ~ ■ 

( = 1 1 !=1 

Finally, to assess the quality of the selection procedure, we compute 
the ratio R,, of the empirical root mean squared error of d„{k„) and 
the minimal empirical root mean squared error of dn{kn)' 



N 1 N 



Rl = Y^ienAKt) - 0f min YidnAkn) - ef. 



lc„e[l,350]'. 
;=1 ' !=1 



Results are presented in Table |2l It appears that R„ is usually "close" 
to 1, except for Weibull distributions. In this case, large values of 
R„ together with large values of jU(fc„) indicate that the optimal fc„ is 
larger than 350 observations. 



5 Real data 

Here, the good performance of the adaptive selection procedure is il- 
lustrated through the analysis of extreme events on a benchmark real 
data set. Nidd river data are widely used in extreme value studies 
||T8l l8|. The raw data consist in 154 exceedances of the level 65 m-^s"^ 
by the river Nidd (Yorkshire, England) during the period 1934-1969 
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(35 years). The N-year return level is the water level which is ex- 
ceeded on average once in N years. Hydrologists need to estimate 
extreme quantiles in order to predict return levels over long peri- 
ods. According to [18J, the Nidd data may reasonably be assumed to 
come from a distribution in the Gumbel maximum domain of attrac- 
tion. This suggests to consider Weibull tail-distributions as a pos- 
sible model for such data. The adaptive selection procedure yields 

kn = 29. The resulting quantile-quantile plot (obtained by plotting 

the points (loglog(n/z),log (X„_,+i^„)) for i = 1, . . . ,k„ - 1) is approx- 
imatively linear (see Figure |Zl), indicating a good fit of the Weibull 
tail-distribution for x > X,j_j We obtained 0„(fc„) - d„{kn) - 0.89. 
One can plug this result in the Weissman-type extreme quantile es- 
timator proposed in |[T2| to obtain 366m^s"^ as an estimation of the 
100-year return level. Note that this result is in accordance with the 
results obtained by profile-likelihood or Bayesian methods, see \Qj or 

m. 



6 Concluding remarks 

In this paper, we introduce a regression model, from which we derive 
a bias-reduced estimator for the Weibull tail-coefficient 6. Its asymp- 
totic normality is established and an adaptive selection procedure for 
kn is proposed. The efficiency of our approach is illustrated in a sim- 
ulation study and on a real data set. However, in many cases of prac- 
tical interest, the problem of estimating a quantile Xp„ = F"^(l - p„), 
with pn < 1 In, is much more important. Such a problem has already 
been studied in [12| where the following Weissman-type estimator 
has been introduced 



^p„(kn) — X; 



n-k,i+l,n 



^ log(l/p.) \'"^'"^ 
Jog(n/fc„)j 



It is, however, desirable to refine Xp^{kn) with the additional infor- 
mation about the slowly varying function i that is provided by the 
LS-estimates for and b. To this aim, condition {Re{b, p)) is used to 
approximate the ratio f "^(1 - pn)IXn-kn+i,nr noting that 

-k„+l,n — F {Un-k„+l,n)/ 
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with Ui^n < ... < Un,n the order statistics of a uniform (0, 1) sample of 
size n, 



Xn-k„+l,n F~^{Un-k„+l,n) 

d {-lo^pnf ^(-logp„) 



(-log(l - Un.k.M,n)Y ^(-l0g(l - Lr„_,„+l,„)) 

d /log(l/p„) 



e ^ / iog(i/p„) \P _ ^ 
fc(log(nA„)) 



Jog(n/fc„)y 



exp 



The last step follows from replacing Uk,^+i,n (resp. En-k„+i,n) by 
(resp. log(n/A:„)). Hence, we arrive at the following estimator for 
extreme quantiles 

/I iM \\^«(''«) r /iog(i/p„)\P _ -■ _ 
^PnC^n) = X„_,c„+i,„ (^f^^^;^ j expp(log(n/fc„)) 

which is similar to the estimator proposed by [|20|1 in the case of 
Pareto-type distributions. Here, the LS-estimators of 6 and h can be 
used after substituting p by the canonical choice -1. The study of the 
asymptotic properties of such an estimator is the aim of [10 1 . 



7 Proofs of our results 
7.1 Preliminary lemmas 

Lemma 3 For all 1 < j <k„ such that kn ^ oo and k„/n — > 0, we have 

T " = 1 + Op( - — / . I uniformly in j. 
log(n/;) \log(nA„)/ 

Proof of Lemma |3l According to Renyi's representation, we have 



n-j+l 



d_ v-" Jn-e-i+ 



fn-e-j+1 

T 
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i.i.d. 



where /y ~ Exp(l). Since 



denoting 



n-j+lr- 



fn- 



e-j+i 



E 



fn- 



e-j+1 



we have, using Kolmogorov's inequality |21 1 (p. 183), that 



, max \Tjn\ > A < 
\i<y</c„ " / 



Var(ri,0 
A2 ' 



A > 0. 



This implies that Ty„ = Op(l) uniformly in /. Taking into account the 
fact that 



" 1 

Y^j- log(n/;) 



0(1) uniformly in;, 1 < j < k„, 



it is easy to deduce Lemma HI 

Let us introduce the £^ -function defined by the integral 



□ 



f ^ , 

:= dt 



for a positive integer m. The asymptotic expansion of this integral is 
given in the following lemma. 

Lemma 4 As x — > oo,for any fixed positive integers m and p, we have 

1 



e~^ \ m m(m + l)...(m + p - 1) 

EUx) = — \i- - + ... + {-ly— — - + ^, , 



The proof of this lemma is straightforward from fT] p. 227-233 and 
the 0-term can be obtained by a Taylor expansion with an integral 
remainder. Denote 

/=1 

The next lemma provides a first order expansion of this Riemman 
sum. 
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Lemma 5 IfK^ oo, kjn ^ 0, j^^^ ^ oo and ^ 0, i/zen 

[ip ~ Cp{lo§{n/k„))~^ as n ^ oo, where Cp = I (log x + < oo. 

Jo 

Proof of Lemma m Denote a„ = j^znjp;- Thenxjt,, can be rewritten as 



^ ZT"^ ZT /«0'/^") - I /nW^^ + I fnix)dx =■— + Ti + T2, 



where = (!-«„ log x)"^ x G [0, 1]. Denoting by fjl\ i e {1, 2}, the 
ith. derivative of /„, we infer that 

Tl = I, {j/kn-t)A'\j/kn)dt 

y, I I {x-t)ff{x)dxdt+ I /„(x)rfx 

J;/lc„ J//lc„ Jo 

=: r3 + r4 + T5. 

Remark that 

ic„-l 



+ 



" \ / ■^'^n Jl/k„ 



2k 



Since /^^' is positive and decreasing on 1 j for n sufficiently large, 
we can prove that 



f^\llK)-ff{l) =o{llK), 
Ts = 0(1A„), 

n\ < =0(1), 

Tv = -^Jfn{l)-fn{l/kn))=0{l/k„), 



13 



and consequently Ti = 0(1 Besides, a direct application of 
Lemma IHprovides 

T2 = l-an + 0{al). 
Therefore X)c„ = 1 - an + 0(1 /?c„) + O(a^). Now, we can check that 



where 

= ^£|aOg(;A.) + 1 + Eny-(\0g{j/k„) + ij 



k„-l 
^" ;=1 



with e„ = O (a„ log^ fc„ j + O (j^^ which tends to by assumption. 

Since Zy=i(log(y/^)5) + 1/ in order to conclude the proof of 

Lemma Is} we only have to remark that — > 0. □ 

7.2 Proof of Lemma[l] 

Remark that 

:= ;log(n/;)[A„(;)-A„(; + l)] 
= 0;log(n/;)log(E„_y+i,„ 

= log(n/;) i{En-j+l,n - E„_j^n)/K-j,n 

where E*_^.jj G [E„_j^„; E„_j+i^„]. Consequently, from LemmalU 

a^,n = efj + Opi- — 7^7-^1 
\log(n/fc„); 

= 0/y + Op(fc(log(nA„))), (10) 

by the assumption x|fc(x)| cxi as x oo with a Op-term which is 
uniform in /. Lemma [His therefore proved. □ 
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7.3 Proof of Lemma m 

We consider 

lij,n := ;log(n/;)[B„(;)-B„(; + l)]. 

In order to study this term, we will use the notations Aiy = En-i+i,nl^n-k„+i,n, 
A2; = En-j,„/En-k„+i,n and yk„ = E„-k„+i,n, and we rewrite jSy,„ as 

^i,n = i log(n/;)|log^(A2; ^ yu,) - log^(A2/y^„)|. 

It is clear that 1 < Aij/Azj — > 1 uniformly in / by Lemma|3]and there- 
fore for n > No, Aij/Aij e [1,2] in probability Under our assumption 
{Re{h, p)) on the slowly varying function, we deduce that 

^j,n = i log(n/;){fc(A2;y,c„)i<:p(Ai;/A2;)(l +Op(l))). 

p 

Now, since A2; — > 1 uniformly in / and h{.) is regularly varying with 
index p, b{A2jyk„) = A'^.b{yk„)i'l + Op(l)) with a Op(l)-term uniform in 
Therefore 

lij,n = i log(n/;)&(y,„){A5^.Kp(Aiy/A2y)(l +Op(l))). 
Again, uniformly in j, 

KpiA,j/A2j) = {A,j/A2j - 1) (1 + Op(l)), 
which implies that jSy „ can be rewritten as follows: 

(}j,n = - j login/ j) biyO{A2j - Aiy)A^:'(l + Op(l)). 
Therefore, we have 

/ log(n/;) \'' 

with a Op(l)-term which is uniform in /. This achieves the proof of 
Lemma 121 □ 
Remark that, since — ^ 1 uniformly in j, one also has 

/ log(n//) 

with a Op(l)-term which is uniform in /, and this proves Corollary [H 
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7.4 Proof of Theorem m 

From model ©, we infer that 

V^(e.ft,)-fl-i(iog(„/fc))i|(ig;5^[) 

^ 1 r- 1 / log(n//) \'' 

= e - - 1) + V^mog(„/)c„)) J- g (jj^) (/,-!) 

(V^fc(log(nA„))). 
Now, an application of Tchebychev's inequality gives that 

Then, under our assumptions. Theorem |2l follows by an application 
of the Central Limit Theorem. □ 



7.5 Proof of Theorem |3] 

From Corollary [T} we have 



log(n/;c„) 
log(n/fc„) fc, 



\log(nA„) " 

Since we have © and ©, the Op-term is negligible. The first term 
can be viewed as a sum of a weighted mean of independent and 
identically distributed variables. Now, using Lyapounov's theorem, 
we only have to show that 



lim — Y EXf = 0, 
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where Xj =[e + fc(log(n/?c„))X;)(l - ^x,„)(/y - 1), ; = 1,...,K and 
s? = L 1 VarX,. We remark that 

Kj, J — L J 

si ~ 0^ V" 1 1 —Xk I as n ^ oo 

^2 "I 

and 

from which we deduce by direct computations that 



9Ca 



by Lemma |5l Our Theorem |3] now follows from the fact that 

si - e%,\og\n/kn). 



□ 



Appendix 

In this appendix, we briefly show how to adapt Hall's class of distri- 
bution function [ , 16,] to the framework of Weibull tail-distributions. 
We introduce the class of distributions D(a, jS) with distribution func- 
tion given by 

1 - F{x) = exp(-H(x)) where H-\x) := + x"^), 

a and j6 being two parameters such that 

< a, < jS < 1 and ajS < 1. (11) 

It is easily seen that under (|TT|) , the above class of distributions fulfill 
assumptions (d) with {Rt(b, p)) and © where 6 - 1/a, p = -jS, i{x) = 
1 -I- x~^ and b{x) = -(ix~^. It is thus possible to obtain distributions 
with arbitrary > and -1 < p < 0. These results are summarized 
in Table m 
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Distribution 


9 


b{x) 


P 


Absolute Gaussian \N\{ii, o^) 


1/2 


1 logx 
4 X 


-1 


Gamma Y{a t 


1 


\osx 

a -a) I 


-1 


WeibulinV(a,A) 


1/a 





— oo 




1/a 







Table 1: Parameters 0, p and the function associated to some distribu- 
tions 



Distribution 





P 


li{k„) 


o{k„) 








K„ 


r(0.25, 1) 


1 


-1 


105.5 


62.2 


1.667 


0.294 


1.26 


186 


r(4,i) 


1 


-1 


222.7 


82.1 


0.548 


0.051 


1.13 


184 


|W|(0,1) 


0.5 


-1 


246.6 


81.1 


0.679 


0.109 


1.21 


189 


nV(0.25,0.25) 


4 


-oo 


305.8 


59.0 


4.016 


0.265 


1.62 


350 


nV{4,4) 


0.25 


-oo 


310.4 


50.9 


0.249 


0.013 


1.43 


350 


2)(1,0.5) 


1 


-0.5 


281.5 


71.1 


0.789 


0.053 


1.14 


43 



Table 2: Simulation results of the adaptive selection procedure 
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Figure 1: Comparison of estimates x x), 0„ (♦♦♦) and 0„ (+ + +) 
the r(0.25, 1) distribution. In (a), the straight line is the true value of 6. 
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Figure 2: Comparison of estimates x x), 6„ (♦♦♦) and 0„ (+ + +) for 
the r(4, 1) distribution. In (a), the straight line is the true value of 6. 
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Figure 3: Comparison of estimates ^4(x x x), G„ (♦♦♦) and 0„ (+ + +) 
the \N\{0, 1) distribution. In (a), the straight line is the true value of 6. 
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Figure 4: Comparison of estimates ^ (with the canonical choice p = -1: 
X X X), e„ (with the true p = -1/2: 000), (♦♦♦) and &„ (+ + +) for the 
D(l,0.5) distribution, hi (a), the straight line is the true value of 6. 
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Figure 5: Comparison of estimates f^6(x x x), 6„ (♦♦♦) and 0„ (+ + +) for 
the 'W{0.25, 0.25) distribution. In (a), the straight line is the true value of 6. 
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Figure 6: Comparison of estimates (^{x x x), &„ (♦♦♦) and 0„ (+ + +) 
the 'W{4r, 4) distribution. In (a), the straight line is the true value of 6. 
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Figure 7: Quantile-quantile plot obtained with kn = 29 on the Nidd river 
data. 
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